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We present a simple method for summing so-called parquet diagrams of fermionic many-body 
systems with competing instabilities using the functional renormalization group. Our method is 
based on partial bosonization of the interaction multi-channel Hubbard-Stratonovich transforma¬ 
tions. A simple truncation of the resulting flow equations, retaining only the frequency-independent 
parts of the two-point and three-point vertices amounts to solving coupled Bethe-Salpeter equations 
for the effective interaction to leading logarithmic order. We apply our method by revisiting the 
X-ray problem and deriving the singular frequency dependence of the X-ray response function and 
the particle-particle susceptibility. Our method is quite general and should be useful in many-body 
problems involving strong fluctuations in several scattering channels. 
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I. INTRODUCTION 

The basic structure of two-particle Green functions of 
interacting quantum many-body systems can exhibit sin¬ 
gularities for different combinations of external momenta 
and frequencies. In simple cases only one particular com¬ 
bination dominates the perturbation series, so that it is 
sufficient to resum the diagrams corresponding to this 
combination. For example, slightly above a supercon¬ 
ducting instability the two-particle Green function be¬ 
comes singular if the total momentum and the total en¬ 
ergy of the two incoming particles are small (particle- 
particle channel). It is then sufficient to approximate the 
effective interaction, appearing in the skeleton diagram 
of the two-particle Green function, by summing only the 
ladder diagrams with small total momentum 1 . Another 
example is the electron gas at high densities, where the 
effective interaction can be calculated by summing the 
geometric series of particle-hole bubbles which dominate 
for small momentum transfer!4 In some cases, however, 
a single dominant scattering channel does not exist and 
singularities in more than one scattering channel appear 
in perturbation theory. In such cases, the usual strat¬ 
egy of summing only particle-particle ladders or particle- 
hole bubbles is inapplicable and one has to solve coupled 
Bethe-Salpeter equations in more than one channel. Di- 
agrammatically, this means summing the so-called par¬ 
quet diagrams where particle-hole bubbles and particle- 
particle ladders are self-consistently inserted into each 
other. Historically, such a resu mmation was first used 
by Sudakov and co-authors 2 ® in the context of high- 
energy meson-meson scattering. Since then, parquet 
methods have played an important role in studying quan¬ 
tum impurity modelJ®®, low-dimensional metalJ^MI, 
liquid HeliunPIEII, and vortex liquid!®. Moreover par¬ 
quet methods have also been used to construct approx¬ 
imations to reduced density matrices of large physical 
system!®, nuclear structure calculation! 20 ^ and differ¬ 
ent lattice models for strongly correlated electrons^DSS. 


In particular, in two seminal papers by Roulet et al . 5 
and Nozieres et alP the threshold exponents of the so- 
called X-ray problenP® were obtained using the parquet 
method. The agreement of the parquet result for the 
threshold exponents with the known exact resullP® at 
weak coupling supports the validity of this technique. 

While parqu et methods are a well established many- 
body too l 31 * 32 ! , they are not straightforward to apply to 
physical problems of interest. One reason could be re¬ 
lated to the fact that the derivation of coupled Bethe- 
Salpeter equations in several channels often relies on 
subtle diagrammatic and combinatoric considerations. 
Moreover, the explicit solution of coupled Bethe-Salpeter 
equations often requires rather substantial algebraic ma¬ 
nipulations or extensive numerical calculations. It is well 
known, however, that the summation of parquet dia¬ 
grams can also be formulated in terms of the renormaliza¬ 
tion group. In fact, the modern formulation of the Wilso¬ 
nian renormalization group for fermionic many-body sys¬ 
tems in terms of a formally exact hierarchy of flow equa¬ 
tions for the one-particle irreducible vertice! 3 ® (we 
shall refer to this approach as the functional renormaliza¬ 
tion group, abbreviated by FRG) opens new possibilities 
for dealing with parquet type of problems. On the one 
hand, the FRG offers a general and systematic frame¬ 
work of deriving parquet equations in a purely algebraic 
way without relying on diagrammatic arguments; on the 
other hand, the formulation of coupled Bethe-Salpeter 
equations as FRG flow equations offers new strategies of 
obtaining approximate solutions. 

In this work we shall revisit the X-ray problem and 
present a systematic approach to obtain the well-known 
parquet result^® within the framework of the FRG. In 
the X-ray problem both the particle-particle and the 
particle-hole channel exhibit logarithmic singularities, so 
that the calculation of the X-ray response function can 
be used as a benchmark for testing truncations of the 
FRG flow equations. Moreover the approximations made 
in Refs. [5] and [B] rely on the logarithmic nature of the 




2 


singularities. However, the necessary algebra is quite de¬ 
manding. Here we show that the threshold exponents 
of the X-ray problem and the frequency dependence of 
the corresponding response function can be obtained in 
a relatively straightforward manner within the FRG. Our 
particular implementation of the FRG relies on the par¬ 
tial bosonization of the theory, in both singular scatter¬ 
ing channels, using multi-channel Hubbard Stratonovich 
transformations. We believe that this technique will also 
be useful in other problems which are dominated by com¬ 
peting singularities in several channels. 

The rest of the paper is organized as follows. In sec¬ 
tion [H] we briefly present the X-ray problem and formu¬ 
late it using functional integrals. In section |Hl| we identify 
the cut-off scheme for the FRG and introduce the regular¬ 
ized deep hole Green function. We write down the exact 
FRG flow equation for the deep hole self-energy and cal¬ 
culate the leading term in the weak coupling expansion of 
the anomalous dimension of the d-electrons. In the suc¬ 
ceeding section [TV] we present our FRG-based approach 
to obtain the particle-hole susceptibility, describing the 
X-ray response, and also derive a similar expression for 
the corresponding particle-particle susceptibility. In the 
concluding section[V]we summarize our obtained results. 


II. FUNCTIONAL INTEGRAL FORMULATION 
OF THE X-RAY PROBLEM 

Before we present the X-ray problem in the functional 
integral framework, we shall briefly describe the nature 
of the problem. The X-ray absorption and emission 
in metals is modeled by the following second quantized 

HamiltoniarPS! 

H = ^2 e k c k c k + e d ( Pd + ^2 U k k'c{c k 'dU, ( 1 ) 

k kk' 

where A k creates a conduction electron with momentum 
k and energy e*. = k 2 /(2m), while d f creates a localized 
deep core electron with energy ed < 0. The Coulomb in¬ 
teraction between the conduction electrons is neglected 
and only the intraband part of the Coulomb interaction 
between the conduction electrons and the deep core elec¬ 
tron is taken into account. For simplicity, the spin de¬ 
grees of freedom are ignored and a separable interaction 
of the form U kk ' = Uu k u k ' is assumed, where the form 
factors u k will be specified below. 

We are interested in the d-electron Green function 

G d (t-t') = -i(T[d{t)d'(t')]). (2) 

Moreover, the experimentally measurable X-ray transi¬ 
tion rate is determined by the Fourier transform of the 
particle-hole response function 

= (3) 

where T is the time ordering symbol and the composite 
particle-hole operator A(t) is defined by 

A(t) = [A(t)d(t) + dHt)c(t)], (4) 


with 


c(t) = ( 5 ) 

k 

The Fourier transforms G d (w) and x ph (w) of the above 
functions for real frequencies to in the vicinity of certain 
threshold frequencies odd and uj x are known to be of the 
fomPUMI 


G d ( W ) 


x ph M 



( 6 ) 

( 7 ) 


where vq is the density of states of conduction electrons at 
the Fermi energy and £ 0 is an ultraviolet cutoff of the or¬ 
der of the width of the conduction band. For small values 
of the dimensionless interaction u = VqU the threshold 
exponents are given by 

? ? = u 2 + G(m 3 ), (8) 

a = 2 u + 0(u 2 ). (9) 

In Refs. [5]and[6] the above results were derived by explic¬ 
itly writing down and solving parquet equations which 
resum the leading logarithmic singularities in both the 
particle-hole and the particle-particle channel. In the fol¬ 
lowing we show that the FRG offers an alternative and, 
in our opinion, simpler way to derive these results. 

In order to set up the FRG for the X-ray problem, it 
is useful to start from an effective Euclidean action de¬ 
pending only on the degrees of freedom corresponding 
to the d-electrons and the relevant linear combination 
([5]) of the conduction electron operators. To derive this, 
consider the Euclidean action associated with the Hamil¬ 
tonian ([I]), 

rP ^ _, 

S= dr^c fe (r)(a T +^fc)c fc (r) 

J ° k 

[P _ 

+ / dTd(r)(d T +£ d )d(T) 

Jo 

fP 

+ U drc(r)c(r)(i(r)fi(r), (10) 

Jo 

where the energies £ k = — A* and = t-d — /r are mea¬ 

sured relative to the chemical potential fji. Although at 
this point we keep the inverse temperature j3 finite, later 
we will take the zero temperature limit /? — > 00 when¬ 
ever it is convenient. The quantities Cfc(r) and d(r) are 
now Grassmann variables depending on imaginary time 
r, and the Grassman variable c(r) represents the linear 
combination c(r) = ^fe u fc c fc( r )- To derive an effective 
action depending only on c(r), we integrate over all c k - 
variables with the exception of the linear combination 
c(t) of interest. Thus we insert the following representa¬ 
tion of unity into the functional integral representation 
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of the correlation functions, 

1=1 V[d, c'] 5 {c\t) - c(t))S(c'(t) - c(r)) 

r 

= J V[c',c'} j V[f],Ti\eio d ^ c '- c )+^-^\ (ll) 


where c' and rj are new Grassmann variables and the 
functional delta-function of the Grassmann variables is 
defined as usuaf 37 L Then we integrate over the c*,- 
variables and subsequently over the auxiliary rj- variables. 
After renaming (/ ->cwe Fourier transform the Gaussian 
part to frequency space and obtain the effective impurity 
model 


Simp = - / (GSM)- 1 ^ - / (Gftiw))- 1 ^ 


rP 

+ U drc(r)c(r)d(r)(i(T), (12) 

Jo 


where the symbol = /3 _1 denotes summation over 
fermionic Matsubara frequencies ico such that for van¬ 
ishing temperature J = f and the Fourier compo¬ 


nents of the fields are defined by dre l “ T c(r) and 

du = Jq dTe lu>T d(r). The Gaussian part of the action 
(12) depends on the non-interacting Green functions 


GSM = £ Sr h ’ (13) 

k 

Go(iu) = —. (14) 

- £, d 


At this point it is convenient to assume that the form 
factor rife is only finite for |£*.| < £o, where £o is a band¬ 
width cutoff of the order of the Fermi energy ep = fi. 
Moreover, we also assume that Uk is such that 


E 


u 


2 

k 


iu: - 



(15) 


where vq is the density of states at the Fermi energy. The 
above integral becomes elementary and we obtain 



FIG. 1. (Color online) (a) G raphical representation of the ex¬ 
act FRG flow equation (191 for the d-electron self-energy. The 
dashed arrows denote the d-electron propagators, the dashed 
arrow with an extra slash denotes the d-electron single-scale 
propagator, the shaded (red) circle with a dot represents the 
scale derivative of the d-electron self-energy, and the shaded 
(blue) square represents the one-particle irreducible effective 
interaction between the d-electrons. (b) Feynman diagrams 
contributing to the effective interaction ( |25[ ) between the d- 
electrons to second order in the bare interaction. The solid 
arrows forming the loops represent the non-interacting con¬ 
duction electron propagators. 


our cutoff scheme. Since the propagator of the c-fermions 
does not exhibit any singularity [see Eq. (16)], we do 
not introduce any cutoff in this sector. On the other 
hand, we regularize the singularity of the non-interacting 
d-electron propagator by replacing 


(G^u ;))- 1 (G^M )- 1 = iux - & + R d A (iux), (17) 

where the regulator function R^iiv) vanishes for A —> 0 
and diverges for A —> oo. The regularized d-electron 
propagator is then 


Ga(*w) 


1 

- £ d ~ £^(*w) + R A {iu) ’ 


(18) 


where the cutoff-dependent d-electron self-energy 
satisfies the exact FRG flow equatiorPSHSU 


d A EiM=/’ Gi{U)T dddd {ux,u'-,J,oj). (19) 

J Uj' 


Gq(jw) = — JTrr'oSgnu; 


1 -arctan 

7r 


(16) 


The effective impurity model defined by Eq. (12) is the 
starting point of our FRG calculation. 


III. DEEP HOLE GREEN FUNCTION 

We begin by discussing the FRG calculation of the d- 
electron (deep hole) Green function. It turns out that the 
leading order term in the weak coupling expansion of the 
corresponding anomalous dimension 77 can be obtained in 
a straightforward way without resumming parquet dia¬ 
grams. To set up the FRG procedure, we need to specify 


Here 


GiM = [-5 a 


( 20 ) 


is the so-called single-scale propagator. The irreducible 
four-point vertex u' 2 \ 0 J 2 , wi) can be identified 

with the properly antisymmetrized effective interaction 
between two d-electrons for a given value of the cutoff 
A. A graphical representation of the flow equation (19) 
is shown in Fig. [l] (a). A convenient regulator for the 
d-electrons is a Litim regulatoi® in frequency space, 


R A (ico) = *sgnw(A — |w|)0(A — |w|), 
so that the cutoff derivative of the regulator is 
d A R A (iuj) = isgnw0(A — |w|). 


( 21 ) 


( 22 ) 
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For simplicity we choose the zero of energy such that the 
d-electron energy vanishes. Ignoring self-energy correc¬ 
tions to the d-electron propagator, we get 


and 


G d K {iuj) 


Ga (iv) 


T-r—— for Iwl < A, 

fcAsgno; 1 1 7 

— for Iwl > A, 

VjJ II 7 


(23) 


isgnu 


e(A-M). 


Since we do not impose any cutoff on the c-fermion prop¬ 
agator, the effective interaction Wj; W 2 , wi) be¬ 

tween the ci-electrons is finite even at the initial scale A = 
A 0 . It should be pointed out that our bare action (10) 
does not contain an interaction of this type. Diagram- 
matically, the effective interaction rQ ddd (o;(,W 2 ,<*d) 
at the inital scale A = A 0 is determined by all diagrams 
with four external d-legs and only c-propagators in inter¬ 
nal loops. In the weak coupling regime, we may approxi¬ 
mate the flowing T'j( ldd uj' 2 ', 0 J 2 , oji) , on the right-hand 
side of Eq. (19), by its initial value at A 0 . For the calcu¬ 


lation of the leading term in the weak coupling expansion 
of the anomalous dimension rj of the d-electrons, we only 
need the effective interaction to second order in the bare 
interaction U. For the special frequencies u)[ = ui 1 = u> 
and W 2 = W 2 = u>' appearing on the right-hand side of 
the FRG flow equation (191 we obtain the effective inter¬ 


action between the d-electrons up to second order in the 
bare interaction 

= -^ 2 [ns c (^-* w ')-n“(o)]+o(t/ 3 ), 

(25) 

where IlQ C (ia;) is the particle-hole bubble with non¬ 
interacting c-fermion propagators, 


ng c (iu)= f G%(iuj - iu)G c 0 (icj). 

J U) 


(26) 


The Feynman diagrams contributing to Eq. (251 are 
shown in Fig. |T] (b). For small |w| we may approximate 
Gq{iuj) ~ — i7rz/ 0 sgnw and obtain 

rj Hwyy.w) = -n{v 0 U) 2 \uj-uj'\+O(U 3 ). (27) 

Substituting this effective interaction into the right-hand 
side of the exact FRG flow equation (19) we obtain a 


closed integro-differential equation for the frequency de¬ 
pendent d-electron self-energy E^(iw). From the numeri¬ 
cal solution of this equation one can obtain the d-electron 
propagator G^(iu}) for all frequencies. Fortunately, the 
singular behavior close to the threshold can be extracted 
from the leading term in the low-energy expansion of the 
self-energy, 


= Sa(0) - (1 - z^)iu + O(co 2 ), 


(28) 


where the cutoff-dependence of the wave-function renor¬ 
malization factor Za defines the flowing anomalous di¬ 


mension. 


(29) 


This quantity can be obtained from the right-hand side 
of the exact FRG flow equation as follows^, 


rjA = ZaA lim 


d 


Oa^a M 


>0 d(iuj) 

= Z K A f G d A {ioj') lim 


d 


0 d(iu>] 


T d A ddd (uj,uj'-,u',u). (30) 


(24) Substituting our perturbative second order result (27) 


and our lowest order approximation (24) for the single¬ 
scale propagator in Eq. p0| we obtain 


(31) 


V = bm ?7 a = (vqU) 2 + 0(U 3 ), 

A—>0 


in agreement with the known weak coupling expansion 
Q. Finally, keeping in mind that our low-energy expan¬ 
sion (281 is only valid for |w| < A, we may estimate the 


d-electron propagator from 


G d (iu>) 


Za=\u\ (M/Aq) 7 ' 


(32) 


Recalling that we have set the threshold energy uj d equal 
to zero, we recover the known threshold behavior (|6j) of 
the d-electron propagator, where we identify £0 = A 0 . 


IV. X-RAY AND PARTICLE-PARTICLE 

RESPONSE 

In this section, using the FRG formalism, we demon¬ 
strate how to obtain the non-perturbative result 0 for 
the particle-hole susceptibility x ph ( w ) which describes 
the X-ray response. Moreover, we shall also derive 
an analogous expression for the corresponding particle- 
particle susceptibility x pp (w). 

To begin with, let us consider the non-interacting 
particle-hole and particle-particle bubbles involving one 
c-fermion and one d-electron propagator, 

II ph M = [ G c 0 {iu)G${iu] - iu), (33) 

J CO 

n pp (iw)= f G c 0 (ioj)G d 0 (-iu + iu). (34) 

J UJ 

Both bubbles exhibit logarithmic singularities. To ex¬ 
tract these, we note that it is sufficient to retain only the 
low-energy part of the local c-fermion Green function, 
Eq. (fl6l) 


Gq(*w) ss -*7Ti/ o sgnw0(^o ~ M). 


(35) 


At vanishing temperature the integrations in Eqs. (33) 


and (34) can then be performed exactly and we obtain 

Co 


IlQ h (ia;) = — v$ In 
IIq P (zcI>) = u 0 In 


iuj + £d 

Co 


(36) 


_iu - 


Va = Ad A lnZ A . 


(37) 
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Due to the divergence in both channels, a simple one- 
channel trunction is not sufficient to calculate the X-ray 
response. The parquet method is a systematic tool to 
resum the leading divergencies in both channels. There 
are different implementations of this method. Our FRG 
formulation is closely related to the particular implemen¬ 
tation of the parquet renormalization group discussed by 
Maiti and Chubukov®, who represent all potentially im¬ 
portant fluctuations via bosonic Hubbard-Stratonovich 
fields, calculate the renormalization of the associated 
three-legged vertices, and finally use these vertices to cal¬ 
culate the susceptibilities. In order to formulate this pro¬ 
gram within the framework of the FRG, we bosonize the 
interactions in the two relevant channels. To do this, it is 
sufficient to retain only those interaction processes which 
mediate the singular interactions in the particle-hole and 
particle-particle channels. To derive the corresponding 
channel decomposition, we note that in frequency space 
the interaction in Eq. (10) can be written in the following 
three equivalent ways, 


ft 

S int = u drc(r)c(r)d(r)d(r) 
Jo 


= U 




/ U)\ J UJ 2 J UJ 2 J Ml 


— U I I I ) (d w / doj '-\-uj ) 

J uj J uj' J Co 

(frequency transfer Cu = lu[ — uji) (38a) 
U I I I ) (d^f C(jj 1 -^-[o ) 

J uo J to' J Co 

(exchange frequency Cu = u>i — uj' 2 ) (38b) 

^ / / / (Qd+cD^— Co) {.d—Uj' Qa/+U>) 

J uj J uj' J Co 

(total frequency ui = ui\ + u 12 ). (38c) 

The last three lines can be generated from each other 
by re-labelling the frequencies. However, if we impose 
an ultraviolet cutoff |cD| < A 0 <C £0 on the bosonic fre¬ 
quency, each of the above expressions describes a dif¬ 
ferent low-energy scattering process: forward scattering 
(38aI, exchange scattering (38b), and Cooper scattering 


(38c). The forward scattering channel can be ignored 


for our purpose because the corresponding susceptibility 
does not exhibit any singularity. We therefore retain only 
the exchange and the Cooper channels for small values 
of the corresponding energies, which amounts to retain¬ 
ing only the following low-energy terms in the effective 
interaction, 

Sint « f 0(A O - M) [~U X A 0 A 0 + U p B a B a ] , (39) 

J uj 

where we have defined the two composite fields, 

Acc — I d^c^j- |_k), (40a) 

J uj 

Bcc = / d—uCu+co • (40b) 


2 Aq 



FIG. 2. (Color online) Low-energy sectors in the three di¬ 
mensional frequency space spanned by the two frequencies uji 
and ui 2 of the incoming fermions and the frequency transfer 
Q = — uii. The blue sector perpendicular to the uj 1 — 012 - 

plane represents the regime |wi + 0 * 2 1 < Ao responsible for 
the dominant renormalization of effective interaction in the 
Cooper channel, while the red sector defined by |u4— Wi| < Ao 
is responsible for the dominant renormalization in the ex¬ 
change scattering channel. 


Although U x = Up = U we have introduced two dif¬ 
ferent coupling constants to distinguish the bare inter¬ 
action U x in the exchange channel from the interaction 
Up in the Cooper (or pairing) channel. Also note that 
a certain small sector of the three-dimensional frequency 
space formed by uii, u >2 and ui = ui’ x — is counted twice 
in Eq. (391, because the sectors defined by \uj' 2 — u>i \ < A 0 
and \ui 2 + wi| < Ao have a finite overlap proportional to 
Aq, corresponding to the intersection of the two frequency 
slices shown in Fig. [2j However, the phase space of these 
processes is proportional to Aq and can be neglected if 
a physical quantity is dominated by generic frequencies. 
Finally, we decouple the interaction using two complex 
bosonic Hubbard-Stratonovich fields x and ifi and obtain 
for the cutoff-dependent effective low-energy action 


Seff 



[(GgM)- 1 ^ + (GgM)" 1 ^] 



\j-J x Xui\ui “b Teo(Au,Xo; 4- A a; ^ a ,)j 
[Up^Qipa + ilpoiB^ipcj + Bo^q)] , (41) 


where the bare values of the Yukawa vertices are 7^0 = 
7 P o = 1 , and we have introduced the notation f A ° = 

/ ffi 0(Ao-M). 

It is now straightforward to write down formally ex¬ 
act FRG flow equations of the above theory for the ir¬ 
reducible verticeP^sj. For our purpose, it is sufficient 
to work with a truncation of these flow equations where 
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only irreducible vertices with two and three external legs 
are retained. Using the same cutoff scheme as in Sec. |III| 
(i.e., we introduce a cutoff only into the d-electron prop¬ 
agator), the cutoff-dependent propagators F*(ioj) and 
{iui) of our bosonic fields are of the form 


II 

o' 

1 

(42) 

Ux 1 + n ph (fw) ’ 

FaM = 

1 

(43) 

Up 1 + n pp (iw) ’ 


where II ph (iu;) and II PP (iio) are the cutoff-dependent irre¬ 
ducible particle-hole and particle-particle bubbles. These 
can be identified with the self-energies of our Hubbard- 
Stratonovich fields. From the formally exact FRG flow 
equation for the generati ng functional of the one-line ir¬ 
reducible verticea^SEfiMHQl ^ j s now straightforward to 
write down formally exact FRG flow equations for the 
vertices of our coupled fermion-boson theory. Neglecting 
the mixed four-point vertices with two bosonic and two 


fermionic external legs (these vertices vanish at the ini¬ 
tial scale and are expected to remain small at least in 
the weak coupling regime) the bosonic self-energies sat¬ 
isfy the truncated FRG flow equations 


9An ph (iw) = / Ga(*w)Ga(*o; — iui) 


x r^ dx (w, ui — ui, w)r^ cx (w — u>, ui, ui), (44a) 
9 aII pp (iw) = — / G c K (iu)G d ^(—iui + iui) 

J CO 

x Y c ^{u), ui — uj,ui)T^(ui — ui,ui, ui). (44b) 


A graphical representation of these flow equations is 
shown in Fig. [3] (a) and (b). The four different Yukawa 
vertices, appearing in Eqs. ( |44a ) and (44b), satisfy the 
flow equations shown graphically in Fig. [3 (c) and (d). 
The flow equations for the Yukawa vertices are explicitly 
given by 


wi,w 2 ) = / -Fa (*w)Ga(-wji + iu})G%{-iuJi +iu) 

J o 


-\cdijj ( / / 




,-AT cd Xl 


d A T d A cx {u } ' 1 ,u 1 ,u; 2 ) = f F^{iu})G A {-iuj[ + iui)G A {-iuJi +*w) 

J CO 

xr^f-wi + w,w(,w)r^(-wi + w,wi, 

SaT^’^ii w 2 ,^i) = - f Fa (iui)G A (iu)[ - iui)G A (iu )2 +*w) 

J CO 

Xr A X ( W 1,^1 - W,w)r^ CX (w2,W2 +w,w)r 

d A r^’(uj i, w 2 ,wi) = - [ F^(iuj)G c K (iu i + iu)G A (iuj 2 - iui) 

J CO 


—+ UI, U} 2 ), 

(45a) 

— w'l + w, 0)2 ), 

(45b) 

W,Wi), 

(45c) 

W, Wi). 

(45d) 


Note that in our cutoff scheme where only the d-electron 
propagator is regularized the FRG equation for the 
fermionic self-energy is still given by Eq. (19). To recover 
the leading order parquet results for the X-ray response 
it is sufficient to set all frequency dependencies of the 
Yukawa vertices equal to zero. In this limit 


tions we find 


<9,II ph (0) = -zz 0 7 1, 
^n pp (0) = I'olp, 

dilx = U)Fp7p7 x , 
da P = -voUx'Yl'yp, 


(47a) 

(47b) 

(47c) 

(47d) 


F 5 dx (0,0,0)=r dcx (0,0,0) = 7„ 
r 5 ^(o, o, o) = r dc ^(o, o, o) = a P . 


(46a) 

(46b) 


where di = —A.d A is the logarithmic scale derivative. Set¬ 
ting u x = VoU x and u p = VqU p and defining g x = 7 ^ and 
g p = 7 p, the last two equations can be written as 

dig x = 2 u p g x g p , (48a) 

dig p = -2 u x g x g p , (48b) 


Ignoring all self-energy corrections (fermionic and 
bosonic) to the right-hand sides of the FRG flow equa- 


which should be solved with the initial conditions g x {l = 
0) = g p (l = 0) = 1. Using the conservation law 
di[u x g x (l) + u p g p (l )] = 0 , we can easily solve the two 
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FIG. 4. (Color online) Graph of the scale depe ndent ver tex 
factors g x {l) = 7^(0 and g P {l) = j p (l) given in (50a 50b I as 
a function of the logarithmic scale factor l = ln(Ao/A) . In the 
inset we show ln^(Z) and In g p (l) which in the leading order 
parquet approximation are replaced by the straight dashed 
lines. 


FIG. 3. (Color online) Diagram (a) represents the FRG 
flow equation ( |44a| ) for the self-energy of the Hubbard- 
Stratonovich field \ (wavy line) representing particle-hole 
fluctuations in the exchange channel. Note that in our cut¬ 
off scheme only the d-electron propagator is regularized. The 
shaded triangles and squares denote the renormalized Yukawa 
vertices, while the other symbols are explained in the caption 
of Fig. [I] Diagram (b) represents the FRG flow equation 
( |44b[ ) for the self-energy of the pairing field ip (zig-zag line). 
In (c) and (d) we show the FRG flow equations (45a 45d[ ) for 
the Yukawa vertices associated with the particle-hole and the 
particle-particle channel. 


coupled differential equations (48a) and (48b) exactly, 


9x( 0 — 


9p{ 0 — 


(u x + Mp)e^“ x+ “ p)! 

(u x -\-Up)l 7 

( u x + Mp)e - ( Ux+ “ p ^ 
u x e( Ua:+u v') 1 + u v e~^ U:c+u p^ 1 


(49a) 

(49b) 


For simplicity, let us now set u x = u p = u, so that the 
flow of the vertices reduces to 


9x(l) 

9p(l) 


e 2ul 

cosh(2 ul) ’ 
e~ 2ul 
cosh(2uZ) 


(50a) 

(50b) 


A graph of these functions as a function of the logarith¬ 
mic flow parameter Z is shown in Fig. [4} Substituting our 
explicit expressions ( 50a 50b ) for the scale dependent 
vertices into Eqs. (47a) and (47b) we obtain 


<9iIP p (0) = -v 0 
diU ph (0) = v 0 


„2 ul 


cosh(2 ul) ’ 

2 ul 


cosh(2 ul)' 


(51a) 

(51b) 


The crucial point is now that our leading order low- 
energy expansion of the frequency dependent susceptibil¬ 
ities is only valid as long as |w| < A = Aoe _i . To obtain 
the leading frequency dependence of the susceptibilities, 
we should therefore integrate the above equations up to 
the logarithmic scale factor Z* = ln(A 0 /|w|). Moreover, 
assuming uZ* < 1 we may repl ace th e fac tor cosh(2wZ) in 
the denominator of Eqs. (51a) and ( 51b| ) by unity. This 
corresponds to the leading order parquet approximation 
discussed by Roulet, Gavoret, and Nozieres 5 . In this ap¬ 
proximation we obtain 


nP h (0) = -v 0 I dle 2ul - "° ri - 2ul * 


vo 

2 u 


n"( 0 ) = 


f 

JO 
1 - 


= ^{l-e 2 
2 u 


Ap 

Iwl 


rl. 


dle~ 2ul = ^°-[l- e - 2 “ r ] 
2u 1 


vo 

2 u 


1 - 


(52) 


(53) 


Finally, we perform an analytic continuation to the real 
frequency axis and obtain for the singular part of the 
retarded particle-hole response function for positive fre¬ 
quencies, 


X P h M = -Re 


nf h (o) 


* 

iuj — 


vo 

2 u 



(54) 


Identifying Ao = £o and keeping in mind that we have set 
the threshold frequencies equal to zero, we reproduce the 
result ([T]) which has first been obtained by MaharPil and 


























































later by Roulet, Gavoret, and Noziered^ using conven¬ 
tional parquet methods. Within our FRG formalism, it 
is straightforward to calculate also the leading singularity 
of the particle-particle response function; we obtain 


X PP M = Re 



iuj— »-ci;+iO+ 


i'O 
2 u 



(55) 


Note that the leading frequency dependence vanishes in 
a non-analytic way for ui —► 0. We have not been able to 
find Eq. (55) anywhere in the published literature on the 
X-ray problem. 


V. CONCLUSIONS 

In this work we have shown how to obtain the threshold 
behavior of the deep hole Green function as well as the 
particle-hole and particle-particle response functions in 
the X-ray problem using the functional renormalization 
group. While our FRG results are equivalent to the lead¬ 
ing order parquet approximatiorP, the calculational effort 
within the framework of the FRG is much lower than in 
the traditional parquet approach. Technically, a novel 
feature of our approach is the use of multi-component 
Hubbard-Stratonovich fields to take into account the sin¬ 
gularities in competing scattering channels. Note that re¬ 
cently several authors have proposed parametrizations of 
the momentum- and frequency-dependent four-point ver¬ 


tex in a purely fermionic formulation of the FRG which 
take the special combinations of momenta and frequen¬ 
cies associated with bosonic collective modes^ 47 into 
account. The method based on multi-channel Hubbard- 
Stratonovich fields proposed here is an alternative to 
these purely fermionic implementations of the FRG. 
In fact, our approach is similar in spirit to the ver¬ 
sion of the parquet renormalization group discussed in 
Ref. [ 2 S], which is also based on the explicit introduc¬ 
tion of fermion-boson vertices via Hubbard-Stratonovich 
transformations. We believe that our method will also be 
useful in the context of other problems where many-body 
interactions lead to competing instabilities in more than 
one channel. For example, in two-dimensional Fermi sys¬ 
tems the low-energy scattering processes of interacting 
fermions can be classified into forward-, exchange- and 
Cooper scattering®^!, which can be bosonized by in¬ 
troducing three different Hubbard-Stratonovich fields. It 
is straightforward to write down the coupled FRG flow 
equations for the irreducible two-point and three-point 
vertices for this problem. The corresponding FRG flow 
will be discussed elsewhere. Finally, let us point out 
that with the FRG it is straightforward to go beyond 
the leading order parquet approximation by taking self¬ 
energy corrections to the propagators and higher order 
irreducible vertices into account. 
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